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SOME EFFECTS OF GEOMAGNETIC ACTIVITY 
ON THE F2 REGION OF THE IONOSPHERE 
OVER LEOPOLDVILLE (4:4°S, 15°3° E) 


By R. G. RAsTOGI 
(High Altitude Observatory, University of Colorado, Boulder, Colorado) 


Received August 16, 1961 
(Communicated by Prof. Sydney Chapman, Hon. F.A.sc.) 


Tue solar diurnal variation of the critical frequency of the F2 layer, foF2, 
has been shown by Rastogi (1959) to be very different at stations having the 
same geomagnetic, but different magnetic latitudes. The value of foF2 
at Ibadan (geomagnetic latitude = 10° N) reaches two maxima in the morning 
and in the evening hours during the years of low solar activity, a feature 
characteristic of the equatorial stations. However, at Leopoldville, the 
diurnal variation of foF2 shows a single afternoon peak, a feature character- 
istic of temperate latitude stations (Rastogi, 1960). These features of foF2 
at Ibadan and Leopoldville are consistent with the magnetic latitudes of 
the stations being 2-5°S and 19°S respectively. 


In a later communication, Rastogi (1961) showed that the lunar semi- 
diurnal variation of midday foF2 at Leopoldville has maximum positive 
devjation at about 11 lunar hours whereas the same at Ibadan is found by 
Brown (1956) to occur at about 03 lunar hours. It was shown that a 
similar type of lunar semi-diurnal variation of foF2 is observed at stations 
having the same magnetic latitudes. 


In view of the large differences in the solar and lunar variations of foF2 
at Ibadan and Leopoldville, a comparison is made in this paper of the effect 


of geomagnetic disturbances in the F2 region of the ionosphere at the two 
stations. 


Martyn (1953) has shown that the disturbance daily variation (difference 
between the diurnal variations on the magnetically disturbed and quiet days) 
of the foF2 and the minimum virtual height, h’F2, of the F2 layer depends 
upon the latitude of the station. There is a general tendency of the electron 
density to be raised at low latitudes but to be decreased at higher latitudes. 


Skinner and Wright (1955) have examined in detail the differences between 
the magnetically quiet and disturbed day variations of the maximum electron 
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density, NmF2 oc (foF2)?, minimum virtual height, h’F2, the height of maxi- 
mum ionization (hmF2)in the F2 layer at Ibadan during the year 1952. 
Following their criterion, the days are taken as magnetically quiet or magneti- 
cally disturbed according to the daily sum of geomagnetic planetary three- 
hour range indices, i.e., 2 Ky < 15 or 2 Ky > 31 respectively. The average 
values of the ionospheric characteristics at Leopoldville for such classified 
groups of magnetically quiet and disturbed days are evaluated separately, 
The virtual height of the F2 layer at a frequency equal to 0-834 times foF?2, 
designated by hpF2, is taken as the height of maximum ionization in the 
F2 layer at Leopoldville. The hpF2 is known to be very close to the hmF2 
values derived by Appleton and Beynon’s method (1940). 


Figure | shows the whole year average diurnal variation of NmF2, h’F2 
and hmF2 or hpF2 on the magnetically quiet and disturbed days at Ibadan 
and Leopoldville. The curves for Ibadan are those given by Skinner and 
Wright (1955). The values of NmF2 at Ibadan are distinctly greater during 
the disturbed than during the quiet days. But there is no significant difference 
between the quiet and disturbed day values of NmF2 at Leopoldville. The 
values of h’F2 for the daylight hours at Ibadan are decreased during the 
disturbed days, but there is only a slight tendency of the same at Leopoldville 
for a few hours around noon. The values of hmF2 at Ibadan are greater on 
quiet than on disturbed days for the forenoon hours, but it is otherwise for 
the afternoon hours. At Leopoldville, there is a tendency for the hpF2 
variation to be opposite to the variation of hmF2 at Ibadan. The disturbed 
day values of hpF2 at Leopoldville are greater than the corresponding quiet 
day values for the forenoon hours, but are less during the afternoon hours. 


The very little difference between the quiet and disturbed day values of 
NmF2 at Leopoldville could be due either to the fact that the magnetic distur- 
bance has little effect on NmF2 at all seasons, or it may have large but opposite 
effects in the different seasons. Therefore, the variations of these parameters 
at Leopoldville are evaluated separately for the northern and southern sol- 
sticial months and are shown in Fig. 2. 


One finds large differences between the quiet and disturbed day values 
of NmF2 during either of the seasons. The electron density, NmF2, for 
any hour of the day is greater on the disturbed than on the quiet days during 
the northern solsticial months (local winter). But, during the southern 
solsticial months (local summer), the values of NmF2 are much smaller on the 
disturbed than on the quiet days. This effect at Leopoldville is similar to 
that at Canberra and Watheroo as found by Martyn (1953), i.e., during the 
~ magnetically disturbed days the mean values of foF2 are depréssed in summer 
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months and are slightly raised in winter months. The disturbance daily 
variation of foF2 at Leopoldville (geomagnetic latitude 6 = 3-1° S, magnetic 
latitude » = 18°S) thus resembles the same at Watheroo (¢ = 41-7°S, 
p= 47°S) and Canberra (¢ = 43-8° S, uw = 47°S), but is very different from 
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that of Ibadan (6 = 10-6° N, » = 2:5°S). This shows the consistency of 
the results grouped on the basis of magnetic latitudes. 


The effect of magnetic storms on the F2 layer can also be described by 
the storm time (Dst) variation of foF2. This is generally expressed as the 


variation of 200 (foF2-foF2)/foF2 with the time since the commencement 
of the storm, where foF2 is the value during the hours following the storm 


and foF2 is the monthly median value corresponding to the same local time. 
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(NmF2), the minimum virtual height (h’F2) and the height of the maximum ionization (hpF2) 
of the F2 layer over Leopoldville during the northern and southern solsticial months of 1952. 
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This quantity approximately equals the percentage variation of the maximum 
electron density, NmF2, in the F2 layer. Matsushita (1959) has studied the 
Dst variation of foF2 at a large number of stations, which were grouped into 
different zones according to the geomagnetic latitude of the stations. It 
was stated that the Dst variations at individual stations within zones so defined 
were very similar to each other. At the equatorial zone within geomagnetic 
latitudes + 10°, the magnetic storm is followed by a short decrease and later 
by an increase reaching a maximum at about 14 days after the commencement 
of the storm. At higher latitudes, the magnetic storm causes an initial short 
increase of foF2 followed by a large decrease with a maximum at about 14 
days after the commencement of the storm. 


In Fig. 3 are shown the storm time variations of foF2 at Ibadan and 
Leopoldville averaged for 19 strong sudden commencement type magnetic 
storms (maximum A, > 50) during the period 1952-1956. The Dst variation 
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Fic. 3. The averaged percentage variation of the maximum electron density: in the F2 
layer at Ibadan and Leopoldville at hours following sudden commencement of magnetic storm 
during 1952-56. 


of foF2 at Ibadan consists of a small decrease during the first day and large 
and consistent increases during the second day of the storm. Contrary to 
this, the Dst variation of foF2 at Leopoldville consists of a short increase 
during the first day followed by large decreases, with the minimum at about 
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40 hours after the commencement of the storm. The curves for the two 
stations are distinctly opposite in nature. The variation at Ibadan is of the 
equatorial type. However, the variation at Leopoldville is of the high latitude 
type. 


It is concluded that the effect of geomagnetic activity on the variations 
of the F2 layer ionization at low latitudes is determined by the magnetic and 
not by the geomagnetic latitude of the station. 
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HETEROCYCLIC COMPOUNDS 
Part XVII.* Synthesis of Dihydrofuro (2': 3'—2 : 3) Quinoline 


By JosEPH P. JOHN AND P. SHANMUGAM 
(Department of Organic Chemistry, University of Madras, Madras-25) 


_ Received August 1, 1961 


(Communicated by Dr: K. N. Menon, F.A.sc.§) 


Furo-(2’: 3’-2:3) quinoline ring system is present in many alkaloids of 
the Rutacee family. Synthesis of dihydrofuro-(2’:3’-2:3) quinolines, 
starting from o-amino ketones and aldehydes has been described in the 
previous paper. The amino ketone or the aldehyde is treated with succinic 
anhydride in hot benzene. The succinoyl derivative, which separates on 
cooling, is then cyclised, by refluxing with 2% aqueous potassium hydroxide. 
The ester of the quinolone acid (II), thus obtained, is reduced to the alcohol 
(II), which is then ring closed to the dihydro-furo-quinoline (IV). This 
paper describes the synthesis of dihydro-furo (2’:3’-2:3) quinoline and 
its 6: 7-methylene-dioxy analogue (VIII). The infra-red spectra of both 
the bases show fairly intense absorption in 1639-1 cm. region. 


Kobashi and Kurvayama‘ obtained the base (IV) by heating 4, 5-dihydro- 
3-furanilide with 85% phosphoric acid at 150-60° for 8 hours. 
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* Part XVI, Proceedings, 1960, 51 A(2). 
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EXPERIMENTAL 
(a) 4’: 5’-dihydrofuro (2':3’-2:3) quinoline (IV) 
2-Quinolone-3-acetic acid (I).—o-Aminobenzaldehyde? (3 g.) and succinic 
anhydride (2-5 g.) were heated together in boiling benzene for 30 minutes. 
The succinoyl derivative that separated on cooling was collected and reflux- 
ed with 2% aqueous potassium hydroxide (250 ml.) for one hour. The 
solution was cooled and acidified with hydrochloric acid when the quinolone 


acid separated. It was collected, washed with water and then with little 
alcohol. It crystallised from acetic acid in needles, m.p.277-—78° C. (dec.). 


Found: C, 64-9; H,4-2. C,,H,NO, requires: C,65-01; H.4-45. 


Ethyl-2 quinolone-3 acetate (II).—The above acid (3 g.) was refluxed 
with absolute ethanol (50 ml.) containing concentrate sulphuric acid (1-5 ml.) 
for 5 hours. The product was collected after dilution with ice-water, and 
crystallised from alcohol in needles, m.p. 186-87°C. 


Found: C, 67:5; H,5-05. C,3H,s;NO; requires: C,67-53; H, 5-63. 


3-(2'-hydroxyethyl)-2-quinolone (I11).—The above ester (1 g.) was reduced 
with lithium aluminium hydride (0-5 g.) in tetrahydrofuran, and the product 
isolated as described previously. It crystallised from benzene in colour- 
less needles, m.p. 196° C. 


Found: C, 69-21; H,5-°7. CyHy,NO, requires: C, 69:8; H5-8. 

4’ : 5’-dihydrofuro (2' : 3'-2 : 3) quinoline (IV).—The above alcohol (0:5 g.) 
was heated with polyphosphoric acid (from 5 g. of P,O; and 3 ml. of phos 
phoric acid) for 4 hours on a steam-bath. After cooling the viscous product 
was poured into ice-water, basified with ammonia and extracted with chloro- 
form, The solvent was evaporated and the residue dissolved in dry benzene 
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and passed through a column of alumina. The benzene was removed and 
the dihydro-furo quinoline was obtained as colourless plates. It crystallised 
from petroleum ether in colourless plates, m.p. 121-22°C. IR. 1639 cm. 


Found: C,76-89; H,5-:10. C,,H,NO requires: C,77-2; H, 5:26. 
(b) 6: 7-Methylenedioxy-4’ :5'-dihydrofuro (2', 3'-2:3) quinoline. 


6: 7-Methylenedioxy-2-quinolone-3-acetic acid (V).—6-Aminopiperonal® 
(5g.) was heated with succinic anhydride (3-5 g.) in boiling benzene for 
30 minutes, cooled and the succinoyl derivative was collected and washed 
with little benzene. It was refluxed with 2% aqueous potassium hydroxide 
(350 ml.) for 1 hour. After cooling it was extracted with benzene to remove 
the free amino-aldehyde. The aqueous layer on acidification furnished the 
quinolone acetic acid. Crystallised from acetic acid, m.p. 290-92°C. (dec.). 


Found: C, 58-33; H,4:1. CygH,NO,; requires: C, 58-3; H, 3-64. 


Ethyl-6: 7-methylenedioxy-2-quinolone-3-acetate (VI).—The above acid 
(3 g.) was suspended in absolute ethanol (100 ml.) and saturated with dry 
hydrogen chloride at 0-5°C. It was kept at room temperature for 24 hours. 
The product, after dilution with water, was collected and crystallised from 
ethanol, m.p. 235-36°C. 


Found: C, 61-08; H,4-73. C,,HisNO, requires: C, 61-09; H, 4-73. 


3-(2'-hydroxyethyl)-6 :7-methylenedioxy-2-quinolone (VII).—The above 
ester (1 g.) was reduced with lithium aluminium hydride (0-5 g.) in tetra- 
hydrofuran and the product isolated as in the previous case. Crystallised 
from alcohol, m.p. 243-44°C. 


Found: C, 61-71; H, 4:65. C,,H,,NO, requires: C, 61-8; H, 4°72. 


6: 7-Methylenedioxy-4’ : 5'-dihydrofuro  (2':3'-2:3)-quinoline (VIII).— 
The alcohol (1 g.) was heated with polyphosphoric acid, on a steam-bath 
for 4 hours. After cooling the viscous reaction mixture was poured into 
ice-water and basified with ammonia. The product separated in fine needles. 
Crystallised from alcohol, m.p. 188-89°C. 


Found: C, 67-61; H,4-30. C,,H requires; C, 66-99; H, 4-18. 
LR. 16397! cm. 
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HETEROCYCLIC COMPOUNDS 
Part XVIII.* Synthetic Experiments in the Quinazoline Series 


By M. KRISHNAN} 
(Department of Organic Chemistry, University of Madras, Madras-25) 
Received August 1, 1961 


(Communicated by Dr. K. N. Menon, F.A.sc.) 


THE attempted synthesis of 2-methyl-4-benzylquinazoline resulted in the 
formation of 2-methyl-4-benzoylquinazoline due to facile oxidation of the 
methylene group.' It would be interesting to study the reactions with sub- 
stituted 4-benzylquinazolines with a view of studying the effect of substi- 
tuents on the rate of oxidation of the 4-benzyl portion of quinazolines. 
However, this would necessitate as starting intermediates with substituents 
in the 3-benzyl position. Such compounds are defficultly accessible. So it 
was decided to study the formation of various substituted 2-benzylquina- 
zolines and to study their behaviour with respect to their susceptibility to 
air oxidation. 

The work presented in this paper relates to the preparation of 4-methyl- 
2-benzylquinazolines with substituents both in the benzene part and 2-benzyl 
part of quinazoline. The following scheme of reactions for the syntheses 
of different substituted 2-benzylquinazolines was based on the method of 
Bishler et 


2-Nitro-3 : 4-dimethoxyacetophenone required for the synthesis of 
4-methyl-2-benzyl-7 : 8-dimethoxyquinazoline di, R=R,=R,=H; 
R, = R, = OCHS) was prepared. 2-Nitroveratroyl chloride was condensed 
with the magnesium ethoxy derivative of diethyl malonate in benzene to 
yield 2-nitroveratroyl diethylmalonate. The oily acyl malonate on refluxing 
with propionic acid for three hours underwent hydrolysis and decarboxylation 
to 2-nitro-3: 4-dimethoxyacetophenone. The reduction of the nitroketone 
to aminoketone was carried out with iron powder and glacial acetic acid. 
The cyclisation of 2-phenylacetamido-3: 4-dimethoxyacetophenone (II, 
R=R,;=R,=H; R,=R,=OCH;) to the corresponding quinazoline 
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was effected by means of ammonia. Similarly, 4-methyl-2-benzyl-6: 7- 
dimethoxyquinazoline (III, R= R,=OCH,; R,=R,;=R,=H) and 
4-methyl-2-benzyl-7-methoxyquinazoline (III, Ry = OCH;; R = R, R, = 
R,=H) were prepared in a similar manner starting from 2-amino- 
4 :5-dimethoxyacetophenone and 2-amino-4-methoxyacetophenone respectively. 
m-Nitrophenylacetyl chloride, o-methoxyphenylacetyl chloride, o-tolylacetyl 
chloride and m-tolylacetyl chloride were condensed with o-aminoacetophenone 
and the phenylacetamides were cyclised to the corresponding 2 : 4-disub- 
stituted quinazolines ae means of ammonia. 


NH, 
cH, 
R 
R Y 
Rs Rs 


The substituted benzylquinazoline was dissolved in alcohol and the 
solution was either exposed to air for several days or air was bubbled through 
the solution for several hours, keeping the volume of the solution constant. 
In all cases the starting materials were recovered unchanged indicating that 
no oxidation took place. 


EXPERIMENTAL 


(1) 4-Methyl-2-benzyl-7 : 8-dimethoxyquinazoline (III, R= R,;= H; 
R, = R, = OCH. 3) 


2-Nitro-3 : 4-dimethoxyacetophenone.—A benzene solution of the mag- 
nesium ethoxy derivative of diethyl malonate (24g.) was prepared from 
3-6g. of Grignard magnesium and 10c.c. of absolute ethanol by standard 
procedure. 2-Nitroveratric acid, prepared according to the method of 
Pschorr and Somuleau,* was converted into the acid chloride. A slightly 
more than equivalent quantity of 2-nitroveratroyl chloride (24 g.) in ben- 
zene (75c.c.) was added to a cooled solution of the prepared magnesium 
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derivative of diethyl malonate during 30 mins. After completion of addi- 
tion, the mixture was refluxed for 1 hour and left to stand overnight. The 
magnesium complex was decomposed with dilute acid and the acid layer 
was. extracted with benzene. The organic layer was dried and solvent 
removed. The crude acyl malonate was obtained as an oil. 


The acyl malonate (35 g.), propionic acid (60g.) and conc. sulphuric 
acid (1-8 g.) were refluxed for three hours. When the evolution of carbon 
dioxide stopped, more of sulphuric acid (20c.c. of 4 N) was added and 
the refluxing was continued for three more hours. The solution was poured 
into water and left overnight. The solid was collected and crystallised 
from alcohol. Yield 15g.; m.p. 91-92°C. Found: C, 53-7; H, 4:83; 
CypH,,NO; requires C, 53-33; H, 4°88. The 2: 4-dinitrophenylhydrazone 
crystallised as orange crystals from glacial acetic acid, m.p. 235°C. Found: 
N, 17:00; C,,H,;0;N; requires N, 17-28. 


2-Amino-3 : 4-dimethoxyacetophenone.—2-Nitro-3 : 4-dimethoxyacetophe- 
none (7 g.) was dissolved in glacial acetic acid (50 c.c.) and iron filings (7-5 g.) 
were added during one hour at the temperature of water-bath. After 
reduction, the amino-ketone was isolated in the usual manner and was 
crystallised from light petrol, when a white crystalline product (5g.) of 
m.p. §3-90°C. was obtained. Found: N, 7:5; CygH,N;,O; requires 
N, 7:2. 

2-Phenylacetamido-3: 4-dimethoxyacetophenone (II, R= R, = Ry =H; 
R, = R, = OCH;).—This was prepared by the reaction of phenylacetyl 
chloride with the above amino-ketone. The amide was crystallised from 
benzene and petroleum ether, m.p. 115°C. Found: N, 4-85; C,,H,,NO, 
requires N, 4:5. 

4-Methyl-2-benzyl-7 : 8-dimethoxyquinazoline (Ill, R= R,=R,=H; 
R, = R, = OCH;).—The foregoing amide (1:5g.) was mixed with dry 
ammonium acetate (15g.) and dry ammonia was passed continuously 
through the molten solution kept at 165° C. for 4 hours. It was then cooled 
and poured into water. The aqueous solution was repeatedly extracted 
with solvent and the extract washed with dilute acid. Basification of the 
acid extract yielded the base. The quinazoline was crystallised from dilute 
alcohol, m.p. 147-48°C. Found: N, C,sHigN.O. requires N, 9-6. 


(2) 4-Methyl-2-benzyl-7-methoxyquinazoline (III, R= R,=R,;=R,=H; 
R, = OCH;) 
2-Nitro-4-methoxyacetophenone.—This was prepared from 2-nitro-4- 

methoxy-benzoyl chloride and diethyl malonate.® Its 2: 4-dinitophenyl- 
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hydrazone formed red crystals from glacial acetic acid, m.p. 225°C. Found: 
C, 47°59; H, 3:0; C,;H,3N,O, requires C, 47:89; H, 3-4. 


(II, R = Rs = Ry = H; 
R,= OCH;).—Treatment of phenylacetyl chloride (3-8 g.) with 2-amino-4- 
methoxyacetophenone® (4-1 g.) in the presence of pyridine yielded 2-phenyl- 
acetamido-4-methoxyacetophenone of m.p. 92-94°C. Recrystallisation from 
alcohol yielded brown granular substance of m.p. 96-97°C. Found: 
C, 72-4; H, 6-1; C,,H,,O,;N requires C, 72-6; H, 6-2. 


4-Methyl-2-benzyl-1-methoxyquinazoline (III, R= = 4H; 
R, = OCH;).—The above amide (2 g.) was cyclised to the base in the same 
manner as described before. The base (1-1 g.) was found to be quite soluble 
in organic solvents. Analytical sample was prepared from dil. alcohol, 
m.p. 100°C. Found: N, 10-4; C,;HigN,O requires N, 10-6. 


(3) 4-Methyl-2-benzyl-6: 7-dimethoxyquinazoline (III, R = R, = OCH;; 
R, = = Ry = H) 


2-Nitro-4: 5-dimethoxyacetophenone.—The above ketone was prepared 
from 6-nitroveratroyl chloride and diethyl malonate in an analogous manner 
as described for 2-nitro-3:4-dimethoxyacetophenone. The yields were 
good when compared to the previous method of preparation. 6-Amino- 
veratrophenone was obtained by the reduction of the above nitro ketone.‘ 


2-Phenylacetamido-4 : 5-dimethoxyacetophenone (II, R = R, = OCH;; 
R,=R,=R,=H).—2-Amino-4: 5-dimethoxyacetonhenone was converted into 
the phenylacetyl derivative. The amide was crystallised from benzene and 
petroleum ether, m.p. 138-39°C. Found: N, 4:8; C,gH,,NO, requires 
N, 4-5. 


4-Methyl-2-benzyl-6 : 7-dimethoxyquinazoline (III, R= R,=OCH;; 
R, = R; =R, = H).—The above amide (3 g.) was cyclised by ammonia. 
The base (2 g.) was crystallised from light petrol, m.p. 140-41°C. Found: 
N, 9°8; CygHigN.O. requires N, 9-52. 


(4) 4-Methyl-2-(m-nitrobenzyl) quinazoline (III, R= R,= R,=R,=H 
R, = N 0.) 


2-(m-Nitrophenylacetamido)-acetophenone (II, R= R, = R,= R,; = H; 
R, = NO,).—m-Nitro-phenylacetyl chloride (5 g.) and o-aminoacetophenone 
were condensed to yield the above amide. Recrystallisation from methanol 
afforded white needles, m.p. 120-21°C. Found: N, 9-26; C,.H,,N.O, 
requires N, 9-39. 
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4-Methyl-2-(m-nitrobenzyl) quinazoline (III, R= R, = R,= R3 = H; 
R, = NO,).—The amide was cyclised to the base by means of ammonia. 
The base was obtained as yellow needles from methyl aclohol, m.p. 132°C. 
Found: N, 14-7; C,gHi3N;03 requires N, 15-0. 


(5) 4-Methyl-2-(o-methylbenzyl)-quinazoline (III, R= R, = R,= R,= H; 
R; = CH;) 


2-(0-Tolylacetamido) acetophenone (II, R= R, = Rg=H; Rg = 
CH;).—o-Methylphenylacetyl chloride was condensed with o-aminoaceto- 
phenone in the presence of pyridine to yield the amide. Yellow crystals 
were obtained on recrystallisation from alcohol, m.p. 69°C. Found: 
N, 5-8; C,,H,,NO requires N, 5-57. 


4-Methyl-2-(o-methylbenzyl) quinazoline (III, R= R, = R, = R,g=H; 
R, = CH;).—The above base was prepared by the cyclisation of the fore- 
going amide. White needles of m.p. 59-60°C. were obtained. The base 
was found to be very soluble in alcohol, chloroform and other organic solvents. 
Found: N, 11:5; C,,HyNe requires N, 11-29. 


(6) 4-Methyl-2-(m-methylbenzyl)-quinazoline (III, R= R, = R, = R,= H; 
R, = CH. 3) i 


2-(m-Tolylacetamido)-acetophenone obtained by the condensation of 
m-methylphenylacetyl chloride and o-aminoacetophenone was straightaway 
used for the cyclisation. The amide was cyclised to the base by procedures 
adopted for the similar cyclisations. The base was obtained as a low melting 
solid and purification proved difficult. So 4-methyl-2-(m-methylbenzyl)- 
quinazoline was converted into picrate and analysed. The picrate was 
crystallised from dilute methanol, m.p. 181°C. Found: N, 14-2; C,3H,,N;O, 
requires N, 14-6. 


(7) 4-Methyl-2-(o-methoxybenzyl)-quinazoline (III, R= R, = R,=H; 
= OCH;) 


2-(0-Methoxyphenylacetamido)-acetophenone (II, R= R, = = Ry= 
H; R,;=OCH,).—The amide was obtained in the usual manner by the 
condensation of 4-5 g. of o-aminoacetophenone and 5-5 g. of o-methoxy- 
phenylacetyl chloride in the presence of pyridine. The amide was crystal- 
lised from methanol; yellow crystals, m.p. 97-98° C., were obtained. Found: 
N, 5:15; C,,H;zNO, requires N, 4:95. 


4-Methyl-2-(o-methoxybenzyl)-quinazoline (III, R=R,\= R,=R,y=H; 
OCH,).—The above amide (3 g.) was cyclised by ammonia in the presence 
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of dry ammonium acetate. Working up the product in the same way, the 
base was obtained as a yellow crystalline substance. It was crystallised 
from light petrol; yellow crystals, m.p. 92-93°C. Found: N, 10-77; 
C,,H,,N,0 requires N, 10-6. 


SUMMARY 


The synthesis of substituted 2-benzylquinazolines is reported. 
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ABSTRACT 


This paper describes an approximate method to calculate the 
incompressible laminar spanwise boundary layer on-an infinite yawed 
uniform cylinder with distributed suction. The calculation for the 
chordwise layer is supposed to be known. A similar method to calculate 
the laminar spanwise layer for a solid boundary is given by Rott and 
Crabtree, but for a porous boundary with suction no such method is 
known to the author at present except the one given in this paper. The 
present method is equally applicable to a solid boundary also. A com- 
parison of the calculation by the present approximate method is made with 
the exact solutions of the spanwise layer with the chordwise velocity 
distribution U = Cx™ at the edge of the layer, and the agreement is found 
to be quite satisfactory. Moreover, the method is very simple in operation. © 


1. INTRODUCTION AND SUMMARY 


THis paper describes an approximate method to calculate the incompressible 
laminar spanwise boundary layer on an infinite yawed uniform cylinder 
with distributed suction. The calculation for the chordwise layer is supposed 
to be known. A similar method to calculate the laminar spanwise layer 
for a solid boundary is given by Rott and Crabtree,’ but for a porous boundary 
with suction no such method is known to the author at present except the 
one given in this paper. The present method is equally applicable to a 
solid boundary also. A comparison of the calculation by the present appro- 
ximate method is made with the exact solutions of the spanwise layer with 
the chordwise velocity distribution U = Cx™ at the edge of the layer, and 
the agreement is found to be quite satisfactory. Moreover, the method 
is very simple in operation. 
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2. NOTATION 


x,y,z = Co-ordinates measured in the tangential direction, spanwise 
direction and normal to the surface respectively (Fig. 1). 


Fic. 1. Sketch showing the orthogonal curvilinear co-ordinates for the flow past an infinite 
yawed, uniform cylinder. 


u,v,w= Velocity components in the boundary layer in x, y and z 
directions respectively. 


Chordwise velocity at the edge of the boundary layer. 
A constant in the expression U = Cx™. 

Index in the relation U = Cx™. 

Spanwise component in the free stream. 


Static pressure in the boundary layer. 


dU 
ax 


Normal suction velocity at the surface. 
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p = Coefficient of viscosity. 


Density of the air. 


Kinematic viscosity. 


displacement thickness for the chord- 
wise layer. 


%£ = f + ( i— v) dz = momentum thickness of the chord- 
wise layer. 


displacement thickness for the span- 
wise layer. 


momentum thickness for the span- 


Vo wise layer. - 
H, = _ = shape factor for the chordwise profile 
ad in the boundary layer. ; 
Hy, = oy" = shape factor for the spanwise profile 
v in the boundary layer. 


mixed momentum thickness of the 
boundary layer. 


dy = a parameter for the family of span- 
wise profiles. 


non-dimensional skin-friction para- 
meter for the spanwise boundary 
layer. 


a shape parameter for the spanwise 
profile. 


ratio of the mixed momentum thick- 
ness to the chordwise momentum 
thickness. 


le, 
af 
co 
* = 1 == = 
v 
* = = 
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I _ By fdv 
0 
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ratio of the momentum thickness of 
the boundary layer in the spa2wise 
and chordwise directions. 


a suction parameter for the chord- 
wise boundary layer. 


a suction parameter for the spanwise 
boundary layer. 


a non-dimensional variable. 


a measure of the chordwise boundary 
layer thickness. 


a measure of the spanwise boundary 
layer thickness. 
a non-dimensional variable in the 


chordwise profile. 


form parameter for the chordwise 
profile. 


a non-dimensional variable in the 
spanwise profile. 


form parameter for the spanwise 
profile. 


ratio of the measure of the spanwise 
boundary layer thickness to that 
of the chordwise. boundary layer 
thickness. 


ratio of the chordwise momentum 
thickness to the measure of the 
chordwise boundary layer thick- 
ness. 


ratio of the spanwise momentum 
thickness to the measure of the 
spanwise boundary layer thickness. 
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ratio of the mixed momentum thick- 
ness to the measure of the span- 
wise boundary layer thickness. 


coefficients in the expressions for 
Sxy- 

a variable to denote Hy; or Hy in 
equation (33). 


a variable to denote K,; or Ky in 
equation (33). 


3. OUTLINE AND DEVELOPMENT OF THE METHOD 


There are three important steps in the construction of this method. The 
first is to recast the mixed momentum equation in a form which is convenient 
for step-by-step computation. The second step is the construction of curves 
or numerical tables for s and q in terms of Hz and Hy, which respectively 
determine the shapes of the chordwise and spanwise velocity profiles. This 
numerical relation between s and q is evaluated by using a family of approxi- 
mate velocity profiles introduced by Schlichting? for the chordwise and span- 
wise layers. The third step is to simplify the calculation of Hy by using 
Iglisch’s solution® on a flat plate with constant suction. The basic assump- 
tion underlying the first and the third steps is that the spanwise profiles are 
approximately a function of suction only on (z/#y) scale. This assumption 
was suggested by the results obtained by the present author.’ 


3.1. Simplification in the mixed momentum equation 


The boundary layer equation for the incompressible laminar flow over 
an infinite yawed uniform cylinder are 


(1) 


(3) 


and the equation of continuity is 


wu, 


fey 
av 
bet = 0,1,2,3 = 
t 
Ky = 
t 
ox dz p ox dz? 
dw 
u—+w— (2 
ox + oz ) 
0 
dz 


176 KAPILDEO PRASAD SINHA 


The boundary conditions are: 


(a) at the surface z= 0, u=v =0 and w= wz is zero or non-zero 
according as the boundary is solid or porous, and (b) at the edge of the 
boundary layer 


z—oo, u=U, v=V, 

and 
dy 
dz oz 


If the flow outside the boundary layer is irrotational 


(5) 
Since the cylinder is infinite in length, 3U/d,=0, so that » = constant 


= V, (say). Also, using equation (3) and the boundary condition at the 
edge of the layer, equation (1) gives 


adU_ 
(6) 


The method really starts from the mixed momentum equation 


{Fu (Vo — v) dz} =y + WoVo (7) 


which is obtained by integrating the spanwise layer equation (2) for a porous 
boundary with suction. Introduce the mixed momentum thickness 


fi (1-y)az (8) 


and the equation (7) reduces to 
d 
(Ulay) = (52), + 0) 


In case of suction there are two arbitrary distributions :—the pressure gradient 
and the suction. As mentioned before in Section 3.0, the spanwise profiles 
on (z/6,) scale will be taken as functions of suction only. Thus a spanwise 
profile is given by 


z 
where ay depends on suction only and consequently so do /y and my. 


4, 

: 
: 
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The ~~ of the velocity profile at the wall is given by 

which transforms equation (9) into 


6, 
(UF xy) ly ry, 


(12) 


_ Wo 


The boundary condition of the spanwise layer equation (2) on the porous 
wall gives 


(13 
) 
which is also a function of suction only. 
Let us define 


6 6. 
= d = 
and then equation (12) becomes 
= (14) 


To find the functional relationship between /, and Ay, Iglishch’s solution® 
for a flat plate with constnat suction is used. The curves of Jy against Ay _ 
is shown in Fig. 2. The values of Jy calculated from the exact solutions 
obtained by the present author‘ for the potential velocity distribution U = Cx™ 
are also shown and they lie approximately on the same curve. The polyno- 
mial given by Thwaites® for this curve is 


= 0-225 + 0-36Ay, + 0°38A,?, (15) 
which transforms the equation (14) into 


(sU0,) = v (0-225 —0 +0- (16) 


This equation for a solid boundary becomes (putting Ay = 0) 
d 0-225 


which is the same as obtained by Rott and Crabtree, 
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ty. 20 SKIN-FRICTION PARAMETER . 


Ay = (SPANWISE SUCTION PARAMETER.) 


XACT 
PLAT PLATE IGLISCH SOLUTION WITH 


WOMOGE NEOUS SUC TION . 
o FLOW NEAR THE STAGNATION LINE 
(61). 
0LUTIONS OF PRESSURE RISE DOMAIN 


or (P<0). 


HLING AND BUSSMANNS 
kESULTS For pro MO We 


> l Ay L 
01 02 O3 


Fic. 2. Variation of J, with yy based on some exact solutions of the spanwise boundary- 
layer with suction. 


It appears that equation (16) is not sufficient to yield the characteristics 
of the spanwise layer, since s, g and Ay are all unknown. Since Ay = gz 
where A, is given from the chordwise layer, only a relationship between 
s and q is required for the solution. 


3.2. A universal relationship between s and q 


In the case of a solid boundary it was analytically shown by Rott and 
Crabtree? that s is a universal function of g and my. This was deduced on 
the basis that chordwise profile can be represented by a function 


mz). (18) 


where mz represents the shape of the profile. The spanwise profile was 
of fixed Blasius profile shape and was given by 


v _ z 
my = 0). (19) 

When suction is applied the shape of the chordwise profile is deter- 
mined by both the pressure gradient and suction and that of the spanwise 
profile is a function of suction only. In the presence of suction, boundary 
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layer profiles are not immediately determined by mz and my, and therefore 
the shape factors Hy and Hy, are chosen to define the chordwise and spanwise 
profiles respectively. Lin* (p.95) also used this parameter for this purpose. 


Thus the chordwise and spanwise profiles can be represented respectively 


The mixed momentum thickness in this case is 


Ss He) [1- Hy) | ae, 


which in terms of s and g and new variable € = (z/#,) becomes 


$= fre H) Hy)] a. 


s(q, Hz, Hy). (21) 


Hence s is a universal function of g, Hz and Hy, and the required relation- 
ship may be evaluated when the family of profiles is specified. 


At present the only system of approximate profiles with suction is that 
due to Schlichting.2 This system of profiles is used here to evaluate the 
relationship between s and q for given Hz and Hy. 


3.3. The calculation of s and q in terms of Hy and Hy 
Thus the chordwise profile is 
1 —e% — sin ne], 


=1-(1+K,)e%, 


where 


> 
8, being a measure of the chordwise boundary layer thickness, 


Nz 


‘ 
and (20) 
0 
‘ 
- 
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The momentum thickness of the chordwise layer i is given in its non-dimen- 
sional form by 


6, 


The displacement thickness is given in the same way by 


The form parameter H, is therefore 


a 1 — 0-09014 K, 
=~ 0-5 +0-06656 Ky — 0-02358 Kz? 


The — velocity profile is given by 
=1—e% + Ky[1—e% — sing ny], 
=1-(1+Kye%, w>3, 


= 0-5 + 0-06656 Kz — 0:02358 Kz? = gp. (23) 


(25) 


_ 2 
5,’ 


5, being a measure of the spanwise boundary layer thickness. 


The momentum thickness of the spanwise layer in its non-dimensional 
form is given by 


Oy 
= + 0-06656 — 0-02358 K,,,?, 
By Ky (27) 


= 


Similarly, the displacement thickness is obtained in its non-dimensional 
form as 


* 
1 — 0-09014 Ky = fy. (28) 
Thus the form parameter Hy for the spanwise profile is 


1 — 0.09014 Ky 
~ +0: 06656 K, — 0: 02358 


— 


(29) 
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and then, the mixed momentum thickness is separately calculated for the 
two cases, y<1 and y>1. 


For y< 1 (i.e., 5y > 5,) the mixed momentum thickness is given by (Cal- 
culation not given here) 


= Sxy 


(30) 


where the coefficients ay, a, dg and a, are given by 


1 


For y> 1 (i.e., 5y> 5,) the mixed momentum thickness is given by 
Dory = by + + b2Ky + = Sry (31) 
where the coefficients bo, b,, b, and b, are analytically expressed as 


2 


1 
a, = 1 — — —__,, 
(8) 
and 
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It is also seen that for y = 1 the two equations (30) and (31) become 
identical. From equations (23), (27), (30) and (31) we get 


— Say 
AY 


(32) 


8y 


If Kz, Ky and y are given, gx, gy and gzy, and consequently s and g can be 
calculated. 


Range of calculation —The functions Kz and Ky determine H, and Hy 
as given in equations (25) and (29) respectively. But for the present cal- 
culations Hz; and Hy are supposed to be given and consequently K, and 
Ky are determined. For this purpose equations (25) and (29) are identical 
in form and may be expressed as equation 


0-02358 HyKy? — Ky [0-09014 + 0-06656 Hy) + (33) 


where K; is a variable representing Kz or Ky. It is known from the results 
obtained by the present author* that the smallest value of H, is 2-0 for the 
asymptotic suction case. This value will also be taken as the lower limit 
of H,, although Thwaites® has obtained some values of Hy less than 2. If 
H; equals the minimum value 2-0, the product of the two roots of Ky in 
(33) becomes zero. One root Kz; = 0, gives the asymptotic suction profile 
and the other root is positive which is excluded by Schliching.? For all 
H; > 2 the product of the two roots of K; is always negative, and only the 
negative roots are permissible. Thus one value of H; yields only one relevant 
value of K;¢. 


The range of Hy is given by 2< Hy < 2-59 where the lower limit is 
fixed by the asymptotic suction profile, and the upper limit by the Blasius 
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profile for the flat plate. It is noted here that for any solid boundary 
Hy = 2°59. The range of H, is taken to be 2<H,;=4. The choice of 
the upper limit is a bit arbitrary, but it is expected that calculations of the 
boundary layer will not be normally taken beyond H, = 4. This limit can 
be extended beyond 4-0 if necessary but it is not proposed to do so here. 
Results of the calculations.—For given values of H, and Hy, Kz and Ky 
are obtained by equation (33). These values of K,, K, together with suitable 
values of y determine s and q with equation (32). Curves of s and q for two 
extreme cases Hy, = 2:0 and Hz = 4-0 for Hy = 2:0, 2:1, 2:2, 2-3, 2-4, 
2-5 and 2-6 are shown in Figures 3 and 4. For given Hy, Hy and s, g can 
be obtained from tables or curves obtained for s and g for given H, and Hy. 


25 


Hx =2:0 


0 as 10 15 


Fie: 3; Curves: for calculating g for H, = 2:0 and arbitrary values of Hy, and s. 
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Fic. 4. Curves for calculating g for H, = 4-0 and arbitrary values of H, and s, 
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3.4. Method for the calculation of Hy 


[t has been mentioned in Section 3.1 that A, depends only on suction 
and may be taken as a suction-parameter. Thus H, depending on suction 
only may be taken as a function of A, only. The relationship between Hy 
and Ay, is now required in the application of the method as obtained in 
Section 3.0. This relationship can be evaluated using Schlichting’s* ap- 
proximate profiles as was done for s and qg in Section 3.3. But it is expec- 
ted to obtain better results by the use of the well-known exact solution 
given by Iglisch* for the flat plate with uniform suction. 


Iglisch* gives a table of A, and Hy, and for convenience, the values are 
reproduced here in Fig. 5. 


260, 


EXACT SOLUTIONS. 


2:50 PAT PLATE IGLISCH RESULTS FoR 
WOMOGENIOUS SUCTION. 


2:40 
DOMAIN 
SCHLICHTING AND BUSSMANN'S 
2:30 RESUTS POR Won 


2:20 


* 


Fic. 5. Variation of the shape-parameter H, with the suction parameter A, of the spanwise 
boundary layer based on some exact solutions. 


3.5. Comparison with exact solutions 


Comparison of q with exact solutions.—In Section 3.3 the universal 
relationship between s and g for given Hz and Hy has been evaluated by 
the use of Schlichting’s? system of approximate velocity profiles. Before 
applying these results to any general calculation on a yawed cylinder, it is 
desirable to estimate their accuracy. For this comparison it is assumed 
here that Hz, Hy and s are given correctly and only the values of g are com- 
pared from the two sources. The typical cases of flow chosen for this purpose 
are:—the flat plate with zero pressure gradient, the flow near chordwise 
stagnation line with constant suction, and some suitable cases of flow with 
adverse pressure gradients and suction for the chordwise velocity distri 
butien U = Cx™(m Z0) at the edge of the boundary layer. 


| 
A FLOW NEAR THE STAGNATION UNE 
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* It is known that on a flat plate either with or without suction, the 
solutions of the two boundary layers (chordwise and spanwise) are identical. 
Hence H, = Hy, s = 1 andg=1. From Fig. 3 it is seen that for H, = Hy 
and s= 1, q is always one. 


For the flow near the chordwise stagnation line with constant suction 
(U = Cx), the comparison is made with the exact solutions and the results 
show that the percentage error is less than 0-5. 


For cases of flow with adverse pressure gradients and suction (U= Cx™, 
m Z0), results of the comparison show that the errors in the calculations 
are less than 6 percent. even when large suctions are applied. 


Comparison of H., with exact solutions.—In this comparison A, is supposed 
given and values of H, from exact solutions are compared with the values 
of Hy obtained from Iglisch.* The same typical cases are compared here 
also. These exact values of H, are plotted against A, in Fig. 5 and are seen 
to lie approximately on the Iglisch*® curve. 


The percentage of error involved is also less than 3 even for adverse 
pressure gradients. 


3.6. Procedure for practical calculation 


To carry out the calculation of the spanwise layer, the results of the 
chordwise layer are assumed to be given. The calculation may start from 
any point, where the characteristics of the spanwise layer are known. Gene- 
rally, it starts from the stagnation point where g and s can be determined. 
If the suction does not start at the stagnation point, Hy is known to be 2-59 
and H, will be already given from the chordwise layer. 


Hence, g and s can be simultaneously obtained by the intersection of 
the curve of s against g for given H, and H, (from the universal relation- 
ship) with the curve sq = 3 [deduced from equation (17) at the stagnation 
point]. If the suction starts from the leading edge, the procedure is slightly 
different. In this case use is made of the exact solution at the stagnation 
point with constant suction. If one uses the exact values of g, s and Hy 
directly as the initial conditions, results are bound to fluctuate for one or 
two steps. This happens so, firstly because the simplified momentum 
equation (16) is not satisfied exactly, and secondly because the characteristics 
of the spanwise layer at this point are not strictly related to the given chord- 
wise results (by an approximate method). Therefore, first of all Hy is fixed 
from the exact solution for given A, which is known from the chordwise 
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flow. Then for given H,; and Hy, g—s curve is drawn according to Section 
3.3. Another curve relating g and s is obtained from (16) near stagnation as 


— 0225 — 0°64 gry + 0-38 
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(34) 


Ay and 


v 


are given from the chordwise flow. The first set of values of gq and s is 
obtained by the intersection of the above two curves. Then Ay = gdyz is 
determined and hence a second value of Hy is obtained from the given 
Iglisch curve. If there is no difference between the two values of 

the initial conditions are satisfactory. If there is any difference, the calculation 
cycle is repeated. In this way values of s, g, A, and Hy, at the starting 
point are determined. 


The calculation then proceeds as follows: The gradient 


at the starting point is calculated by equation (16). The quantity sU0,, at 
the next point x, is obtained from the relation, 


= (SUP x) pax, + (X1 — Xo) | (35) 


Thus s is known at the point x,, since U and 0, are given. To determine 
q at x, Hy is tentatively taken to be the same as at x9. This procedure is 
supported by the fact that H, is constant near stagnation (exact solution). 
q is then obtained from the results of Section 3.3 for known Hg, Hy and 
s. Ay is calculated by the relation A, =A, which again determines Hy 
from the Iglisch* table or curve. If there is any difference between the two 
values of Hy, the cycle of calculation is repeated, and if not, the calculation 
proceeds to the next step. 


When more than two values of H, become known, tentative value of 
Hy, may be chosen from graphical extrapolation. For a solid boundary 
H, is constant (Hy = 2-59) and no extrapolation is necessary. 


It has been found in practice that one cycle of calculation is usually 
enough, and in no case were more than two cycles required. 
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3.7. Discussion and conclusion 


The present approximate method for calculating the laminar spanwise 
layer on an infinite uniform yawed cylinder is mainly constructed for a porous 
boundary, but it may also be applicable for a solid boundary. The method 
has given very good agreement with the exact solutions of the spanwise 
layer for the chordwise velocity distribution U = Cx™ at the edge of the 
layer. The results of the chordwise layer are assumed to be known. This 
assumption of exact knowledge of the chordwise layer is an ideal which 
will not be realised in practice. Any error involved in the approximate 
calculation of the chordwise layer will be automatically inherited in the 
calculation of the spanwise layer. In addition to this unknown source of 
error, the method may be expected to be fairly inaccurate after Schlichting’s? 
separation at 


_ _ 9.0682. 
Vv 


This inaccuracy is the intrinsic weakness of Schlichting’s? profiles for 
which there appears no remedy atthe moment. However, no approximate 
method at present gives very satisfactory results in regions of increasing 
pressure. But it should be noted that the calculation may be continued 
even after this so-called separation point, if necessary. 
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1. INTRODUCTION 


A THEORETICAL analysis is made of the flow of an incompressible non-New- 
tonian viscous fluid contained between two torsionally oscillating infinite 
parallel discs on lines similar to those of the corresponding Newtonian case 
considered by Rosenblat. We study the two cases where (i) one disc only 
is oscillating torsionally while the other is at rest and (ii) both the discs 
oscillate with the same frequency and amplitude but the motions at any 
instant are in opposite directions. Assuming that the amplitude of the 
oscillation namely (2/n) is small and neglecting the terms of 0 (2/n)? in the 
equations of motion in comparison with terms of 0(2/n) we solve exactly 
for the velocity field and pressure satisfying all the boundary conditions 
completely. This approximation gives the same transverse velocity field as 
in the Newtonian case obtained by Rosenblat. We find that the radial- 
axial flow has a mean steady component and a fluctuating component of 
frequency twice that of the oscillating disc. We observe that the approxi- 
mations made gives the velocity field valid only for small Reynolds numbers 
unlike the Newtonian case, where the asymptotic expansion of the velocity 
field for large Reynolds numbers has also been obtained. We have studied 
the non-Newtonian effects on the streamlines of the steady part of the 
radial-axial flow and on the steady component of the radial velocity in both 
the cases through a dimensionless parameter S = (v¢/d*), where v¢ is the 
kinematic coefficient of cross-viscosity assumed as constant and d is the 
distance between the plates. We have compared the flow of the Newtonian 
case (S = 0) with the non-Newtonian case considering both negative and 
‘positive values of S. We observe the continuity of the effect as S changes 
from negative to positive values. We have derived the expressions for the 
shearing stresses acting on the plates for small Reynolds numbers. 


1. We consider a body of non-Newtonian fluid bounded by two infinite 
parallel plane discs which are represented by the planes z= 0 and z=d 
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in a cylindrical polar co-ordinate system. The discs perform torsional 
oscillations about the axis, r=0. If w, v, w are respectively the radial, 
transverse and axial velocity components, p, the pressure, p, the density, 
Jt, He, the coefficients of viscosity and cross-viscesity and v, vg, the kinematic 
coefficients of viscosity, the equations of motion of a non-Newtonian fluid 
in cylindrical polar co-ordinate system are: 


ve { dr dr? + 
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ow ow ow 


Ww ou 1 du du 
du dw + ~) } 
~ dz? ' drdz 


« 


2. ONE Disc OSCILLATING 


We consider the disc at z=0 to perform torsional oscillations of 
frequency » and angular speed 2, while the disc at z = d, remains at rest. 


Following Rosenblat, we assume the velocity field as 


v =rQe'tg(y), 


Further we assume that 


The equation of continuity (1. 4) is identically satisfied by such a choice 
of velocity field. The boundary conditions | of the problem are 


uxw=0, = Re(rQeinty z=0, (2.3) 


190 
{ and 
uu, 
| | (1.4) 
u ay FO: 7), 
(2.1) 
w =— | 
where 
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and 
u=v=w=0 at z=d. (2.4) 
Substituting (2.1) and (2.2) in the equations of motion nai 1) to (1.3) 
we have 
i 
+(2 ) (=) ]- le 
3) dey? 
Q\? dg 
2(* [ 
ig(y) + (5) 
1 Q\? .fdg dF 
of 
st. 4(2)" oF 
1) 
2 °F dK Q\? dF 
2) 
ice 


oF 
F= g=0 aty= Il. (2.9) 
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On the assumption that the amplitude of the oscillations namely Qjn 
is small, retaining only the first order terms in 2/n, we have for the transverse 
velocity component 


(2,10) 


with 
g(0)= 1 and g(l) = O. 
The solution of (2.10) is 


_ sinh ViR (1 — y) 


sinh ViR 


g(y) 


same as obtained by Rosenblat in the Newtonian case. 


For small Reynolds numbers we have 


—y) {ft + 8y — 12y? + 3y%)| cos nt 


+B y—y)sin nt} + 


Noglecting the terms of order (2/n)?, we have from (2.5) for the radial-axial 
component of velocity 


2 3 2 
— = — P(y, + [er (2.12) 


From (2.7) we have 


* (2.13) 


and equating the coefficient of r? to zero we have 


12P sip de 


which determines P(y, 7). 


Let 


F (y, 7) = f(y) + e*, 


| 
0, | (2.14) 
(2.15) 
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and 
K (y, 7) = Ky + Ki (0) (2. 16) 


From (2.14) we have 
P ir _ 


Substituting (2.15), (2.16) and (2.17) in (2.12) and (2.13) we have 


2ih’ (y) — 


= — 38[er + 2.18) 


2ih (y) et” LF" +H" (9) + 


+ K, (y) e%, (2.19) 


These equations are sufficient to determine the unknowns f(y), h(y), K, () 
and K,(y). The above equations indicate that the radial-axial flow has a 
mean steady component and a fluctuating component of frequency twice 
that of the oscillating plate. From (2.11) we have 


[e = E Ay, — cos it Ad +i)y—1 ] 


cosh A — cos A 2 L coshaA(i +i) —1 


(2.20) 


4 cosh A — cos A 


cosh A(i +i) —1 


Mn 
+ 
_| 
and 
| 
and ' 
where 
: 
15) A=VY2R and 
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Substitution of (2.20) and (2.21) into (2.18) yields 


= by + 383 [ 


cosh Ay, + cos 1 [= Ay,—cos Ay, 
~ cosh A — cos A cosh A — cos A } 


(2.22) 


l 


0) — (y) 


=1, +38" 


cosh + i) — 1 


~ 2 | cosha(l +i) — 1 


where we have taken 


2a(r) = Ly + 


The boundary conditions (2.8), (2.9) reduce to 
f=f'=0 aty=0 and =y=1, (2.24) 


h=h'=0 aty=0 and y=1. (2.25) 


The solution of (2.22) for the mean steady component under the con- 
ditions (2.24) is found to be 


$y) = = (sinh Ay, + sin Ay,) 


+ wal — y)*(cosh:A + cos — 3y2 + 2y3) 


x (sinh A+ sin A) — 5 


x [(2y — 3y? + 1) (sinh A — sin A) — (sinh Ay, — sin dy.)I} 


r2 
— 38S — (2. 26) 


and f 
and 
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and 


1° 0) = {[ 4 Cost + 005 ay.) 


+A — y) (sinh A + sin +301 — 4y + 3y?) 


xX (cosh A + cos A) -- 572 3y) 


A2 6 
+ [ cosh Ay, — cos + 1) 


(sinh A — sin + 38 (4y—1—3y). (2.27) 


For small Reynolds numbers we have 


FY) — — 9) + [- Sy> + 35y4 


— 105y* + 175y? — 133y — 21] 


+ BS — syt + toy — 10)]} + 0(R5), (2.28) 


and 
£0) — {6 — + 


[42—1365y2+45 (—y* + 


BS — 5y* + — 10y) + 20}} + 0(R9). 
(2.29) 


The solutions for S = 0, corresponding to the Newtonian case are in 
agreement with those of Rosenblat. 


The unsteady component of the radial-axial flow is given by the equation 
(2.23). We have 


4 [cosh A (1 +i) — 1} feosh 1-49 


) 
) 
4 
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+ ‘) [va + cosh A(1 + 


sinh ( (cosh +) 1) 


ACL +i), 


sinh A(1 + 7) cosh V2 


j) Sin A (1 + +f1 cosh Dy] 


A(l +i) l 


sinh A(1 + i) 
ACL +3) 


1 + cosh + 7) cosh 
A(l 
/2 


x sinh A(1 + 7) sinh + 


x cosh 


x (1 +) — 


sinh A(1 + i) sinh {y— i) 


x \- sinh A(1 + i) y, + sinh A(1 + i) cosh 
/2 cosh A(1 + i) sinh yt 


For small Reynolds numbers we get, 


h(y) ~ y) — 3iRS) 


— 150y? + 167y + 29)} +0(R%). (2.31) 
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The shearing stresses for the velocity field (2.1) are given by 


du dv 
dz 


Tre = fhe 


For small Reynolds numbers the shearing stresses acting on the plate 
z=0 are given by 


2 
+ cos 2nt + R(3S + 5) sin 2nt | + 0(R%, 


Te = + 08 nt nt] +0(R%), 


= — be 40 [ cos 3nt + 3 cos nt 


+R (38 — (sin 3nt + sin nt) | + 0(R%). 


The shearing stresses acting on the disc z= d are 


Tre [1 + cos 2m R(38 + sin 2nt | + 0(R9). 


0 
Tye = — 452 cost + sin nt] + 0(R, 


and 


2 
Trp = — Pe = [cos 3nt + 3 cos nt 


+R (3s + a0) (sin 3nt + sin nt) | + 0(R%). 


du u 
Tre = => | — 2pe = 
oz 
ov ou 
> 
and 
| 
q 
a 
= 
- 
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Figure | represents the effect of non-Newtonian term on the typical 
streamlines of the steady component of the flow. We have drawn the curves 
for S = — 0-5, — 0-1, 0 and +0-05. We observe the continuity of the 
effect as the parameter ‘S’ changes from negative to positive values. The 
negative non-Newtonian coefficient flattens the streamline while the positive 
coefficient makes it more curved. Figure 2 shows the dimensionless steady 


Fic. 1. One disc oscillating: Streamlines of steady radial-axia! ficw fer R=5 and 
S=—0-°5, —0:1, 0, +0°05, 


O06 -O4 -02 
Fig. 2. One disc oscillating: Steady radial velocity f’ for S= — 0-1, 0, + 0-05, 
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component of the radial velocity, namely, f’ for S = — 0-1, 0 and + 0-05. 
We find f’ = 0 at y = 0-45 for S = — 0-1 and at y = 0-48 for S = + 0-05. 


3. Discs OSCILLATING 


We now consider the case where the two discs are oscillating with same 
frequency and angular speed but in opposite directions. The boundary 
conditions in the present case are 

u=w=0, v= Re(rQeinty at z=0, 


and (3.1) 
u=w=0, v=—Re(rQe™) atz=d. 


The equations of motion (1.1) to (1.3) can again be transformed by the 
choice of the velocity field (2.1) in a similar way and for the transverse flow 


we have 


2 
(3.2) 
as previously but with g(0)=1, g(Ij)=—1. 


The solution now is 
_ sinh JiR y, — sinh ViR y 
sy) = sinh 


(3.3) 


once again the same as that obtained by Rosenblat. 
For small Reynolds numbers this yields the velocity profile 


+Bya — y)sin ne} + 0(R%). (3.4) 
From (3.3) we have 


= [cosh (1 — 29) — 0055 (1 2) 
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[cosh (1 + i)(1 — 2y)— 
| (3.5) 
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A A 
cosh + cos 5 
2 


A A 
cosh A — cos A [cosh (I 2y) -- cos > (1 2y)| 


cosh 4 (1 +i)+1 
2 cosh +i) — 1 


x [cosh (1 +00 1] ew. (3.6) 


Substitution of (3.5) and (3.6) in equation (2.18) yields the cquations 
necessary for determining the steady and the fluctuating components of the 
radial-axial flow under the same boundary conditions (2.24) and (2.25). 
The steady components of the axial and radial flows are 


A A 
cosh 5 COS 5 


A 
x (5 sinh 5 — cosh 5) + sinh 5 2y) + y cosh 5 


~ sinh 3] + (384+ 1) — 298) (cos 4 — sin *) 


x Sin 5 (I — 2y) COS 5 +5 sin 5]} (3.7) 


A A 
1 cosh 5 + cos 5 
=3 {(383 1) [6 ») (5 sinh 5 


— cosh >) cosh + cosh (1 +(3s4 + 


x [oo —y) (cos — sin 3) cos 
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For small Reynolds numbers we have 
R R? 


— + 2y — 1)] — tay? — ay + 3)} + 
(3.9) 


+ 5y(13y — 1) ~ — BS cayt— 


+ Ty? 3y) + 3] + 0(RS). 


The fluctuating component of the radial-axial flow is given by 


+4) 


V2 h(y) 


A? i 
38 + 5 
coshA (1+ i) —1 


x y + sinh XC + i) cosh 3 (1 +i) 


Ad +i, 


sinh2 
sinh sinh V2 


AU +i) 


(383 - 5) 


+i) ~ sinh 51+) -2y)} , 


(3.11) 
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[cosh 1] [2 cosh “ 


A(1 + i) A(1 +i) 
sinh 


osh 4 +i) sinh 5, a+at 


+ (38% + i) cosh 


x +i)+ i) (1 osh 
+) 


D [cosh + i) — 1} {2cosh i) 


+i) 


feosh + i) ¢ — cosh 


+ sink 5 + ome: + - +5) 


x {eosh?} (1 + i) (1 — cosh on MED i 


+ sinh 5 ( +i) cosh 5 sinh +ah. 


(3.13) 


In the present investigation we find that our approximations are valid for 
small Reynolds numbers only... Figure 3 shows the_effect of non-Newtonian 
term on the typical streamline of the steady flow when both the discs 
are oscillating in the opposite directions. The flow is symmetrical about 
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the plane y= 4. Curves are drawn for R = 10 and S = — 0-5, — 0:1, 0, 
+0-05. The effect of negative non-Newtonian coefficient is to flatten the 
streamline. and that of the positive non-Newtonian coefficient is to make it 
more curved. Figure 4 represents the steady component of the radial flow 


18 - $2-05 
282- $=-0-1 
3 83- $2 

4 84’ - $2+0:°05 


Fyg. 3. Two discs oscillating: Streamlines of steady radial-axial flow for R = 10 and 
S =—0-5, —0-1, 0, + 0-05. 


__, FG. 4. Two discs oscillating: Steady radial velocity f fcr S = — 0-1, 0, + 0-05, 
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for $= — 0-1, 0 and +0:05. For S = — 0-1 it is zero on the planes 
y = 0-26 and y=0-72. There is a radial outflow between these and the 
walls due to the predominance of centrifugal action. For 0-:26< y< 0-9 
there is an inflow as the pressure gradient is then the main factor. Foy 
= + 0:05 itis zero on the plane y = 0-29 and y= 0-71. 


The shearing stresses acting on the discs z= 0 and z=d are equal 
in magnitude but opposite in sign. 


For small Reynolds numbers we have the shearing stresses, 


Tre = [4 + R (128 — 1) sin + 0(R4), 


R 3 
Toyz = — 52 [(: + nt — sin nt | + 0(R5), 


r? Qs 


Tre = — He [s cos nt + R (128 — sin nt 


+ R(12S — 1) sin 3nt] + 0(R%). 


SUMMARY 


The flow of an incompressible non-Newtonian viscous fluid: contained 
between two torsionally oscillating infinite parallel discs is investigated) 
The two specific cases studied are (i) one disc only oscillates while the other 
is at rest and (ii) both discs oscillate with the same frequency and amplitudé 
but in opposite directions. Assuming that the amplitude of oscillation, 
Qjn, is small and neglecting the squares and higher powers of 2/n, the 
equations of motion have been solved exactly for velocity and pressuré 
satisfying all the boundary conditions. The effect of both positive and 


flow has been represented graphically. 
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